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1.  Long-Term  Goals 

In  terms  of  target  detection  and  classification,  scattering  from  exposed  rock  on  the  seafloor, 
(i.e.,  individual  rocks  and  rock  outcrops)  presents  some  of  the  most  difficult  challenges  for 
modern  MOM  and  ASW  sonar  systems  in  shallow  water.  Work  on  characterizing,  model¬ 
ing  and  simulating  mean  levels,  and  other  statistical  measures  of  acoustic  scattering  from 
rocks  and  rock  outcrops  is  therefore  critical.  Unfortunately  (and  curiously)  information  on 
scattering  from  underwater  rock  and  outcrops  is  almost  non-existent.  Scattering  from  rock 
outcrops  is  not  simple  enough  to  be  encompassed  by  a  single  scattering  strength  curve,  but 
has  a  variety  of  expressions  depending  on  the  exact  geomorphology  of  the  rock.  Smoothed 
surfaces  may  actually  scatter  less  than  surrounding  sediment;  curvature  may  dramatically  af¬ 
fect  scattering  and  rough  areas  as  seen  on  the  rock  outcrop  in  Fig.  1,  display  high  variability 
which  could  pose  difficulty  for  target  detection  and  classihcation  systems. 

The  primary  long-term  goal  of  this  research  project  is  to  increase  understanding  and 
modeling  capabilities  for  high-frequency  acoustic  scattering  from  rock  and  rock  outcrops.  In 
addition  to  an  increase  in  basic  understanding  of  the  characteristics  of  scattering  from  rock, 
any  resulting  advances  in  modeling  would  be  useful  for  improving  simulation  capabilities  and 
for  improving  detection  and  classihcation  tools.  Inverse  models  based  on  forward  models 
would  be  essential  for  using  sonar  systems  for  remote  sensing  of  seahoor  properties.  An 
understanding  of  spatial  coherence  functions  for  isotropic  and  anisotropic  rough  seahoor 
surfaces  could  allow  a  method  for  separating  natural  ‘target’-like  objects  such  as  rough  rock 
from  man-made  targets. 


2.  Objectives 

Our  objectives  for  the  proposed  study  of  scattering  from  rocky  seahoors  and  rock  outcrops 
are  intended  to  address  many  of  the  open  questions  which  exist  for  scattering  from  these 
types  of  surfaces  and  include  increasing  our  basic  understanding  of: 

(1)  geoacoustic  characteristics  of  rock  relevant  to  scattering, 

(2)  scattering  strength  versus  grazing  angle,  and 

(3)  scattering  statistics  of  rock  outcrops. 
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Figure  1.  Photo  of  rock  outcrop  in  the  Oslofjord  near  Larvik,  Norway.  The 
outcrop  is  similar  to  underwater  outcrops  from  which  acoustic  scattering  was 
collected. 

These  goals  were  achieved  through  examination  of  existing  literature,  analysis  of  held  data 
or  lab  measurements  and  by  the  use  of  extended  approximate  or  numerical  scattering  models. 

3.  Work  Completed 

3.1.  Geoacoustic  Properties  of  rock  outcrops.  As  virtually  no  information  exists  on 
scattering  from  rock  outcrops,  we  have  worked  on  obtaining  relevant  physical  characteristics 
of  rock  outcrops,  such  as  the  geoacoustic  properties,  roughness,  and  morphology  for  use 
in  models  of  acoustic  scattering  from  rock.  This  year  the  methods  developed  in  previous 
years  to  estimate  scattering  strength  from  SAS  data  were  applied  to  a  wider  set  of  data  and 
reported  in  [1]. 

The  historical  literature  was  examined  to  determine  the  mineral  composition  of  the  bedrock 
in  the  area.  These  data  [2]  were  used  to  estimate  the  compressional  wave  speed,  shear  wave 
speed,  and  density  of  rock  outcrops  based  on  Hashin-Walpole-Shtrikman  bounds  [3].  These 
measurements  in  conjunction  with  effective  medium  theory  resulted  in  estimates  6393  m/s 
for  the  compressional  wave  speed,  3276  m/s  for  the  shear  wave  speed,  and  2708  kg/m3  for  the 
bulk  density.  To  provide  estimates  of  interface  roughness,  an  experiment  was  performed  in 
May  2013  using  stereo  photogrammetry.  Both  of  these  types  of  environmental  ground  truth 
were  used  to  compute  approximate  scattering  models  that  are  compared  with  measured  data. 
Analysis  of  these  data  reported  in  [1]. 

Roughness  results  from  the  May  2013  experiment  are  presented  in  Fig.  2  and  Fig.  3  of 
abraded  and  plucked  surfaces  respectively.  The  terms  abraded  and  plucked  are  discussed 
below.  The  root-mean-square  (RMS)  roughness  of  the  abraded  surface  is  2.09  mm,  and  45 
mm  for  the  plucked  surface.  These  quantities  correspond  to  values  of  0.86  and  18  respectively. 
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Figure  2.  (color  online)  Rongh  interface  results  from  a  glacially  abraded 
surface  in  (a)  the  low-resolution  mode,  and  (b)  the  high-resolution  mode.  The 
glaciers  flowed  in  the  negative  y  direction.  The  color  bar  corresponds  to  height 
reference  to  the  surface  mean,  and  the  brightness,  or  black/white  information 
communicates  the  surface  slope.  The  dashed  box  in  (a)  indicates  the  portion 
of  the  surface  that  is  shown  in  (b). 

More  information  on  the  measurements,  as  well  as  parameters  extracted  from  the  power 
spectra  can  be  found  in  [1]. 

3.2.  Scattering  Strength  estimates  from  rock  outcrops.  The  sonar  data  analyzed  was 
collected  in  April,  2011  during  a  joint  held  experiment  that  took  place  near  Larvik,  Norway, 
as  part  of  a  collaborative  work  with  the  Norwegian  Defence  Research  Establishment  (FFI). 
The  SAS  system  operated  at  a  center  frequency  of  100  kHz,  has  a  bandwidth  of  30  kHz  and 
was  operated  from  the  HUGIN  Autonomous  Underwater  Vehicle  (AUV). 

SAS  images  of  a  smooth,  hat  rock  outcrop  can  be  seen  in  Fig.  4  and  is  denoted  R, 
and  SAS  images  of  four  other  rock  outcrops  can  be  seen  in  Fig.  5.  These  outcrops,  called 
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Figure  3.  (color  online)  Roughness  measurement  results  for  the  plucked 
rough  interface.  The  colorscale  communicates  surface  height,  and  the  sur¬ 
face  slope  has  been  included  as  changes  to  grayscale  value  to  accentuate  low- 
amplitude  roughness  not  resolved  by  the  color  scale.  Only  the  low  resolution 
measurement  is  presented  because  the  high  resolution  mode  does  not  include 
enough  facets  to  obtain  a  proper  sample  size. 


roches  moutonees,  were  formed  through  glacial  erosion  and  exhibit  two  contrasting  roughness 
characteristics.  The  stoss  side  is  formed  through  the  mechanism  of  glacial  abrasion  [4]  where 
the  glacier  flows  up  onto  the  outcrop  and  is  characterized  by  a  very  smooth,  polished  surface 
at  small  scales,  and  a  gently  undulating  surface  at  large  scales.  The  leeward  side  is  formed 
through  the  mechanism  of  glacial  quarrying  [5]  where  the  glacier  flows  off  of  the  outcrop  and 
is  characterized  by  a  stepped  appearance,  where  the  step  sizes  and  orientations  result  from 
the  internal  joint  organization  of  the  bedrock,  at  scales  of  0(lm)  and  small-scale  isotropic 
roughness  resulting  from  fracturing  at  small  scales  (O  <  0.1m).  Boxes  denote  areas  from 
which  scattering  strength  was  estimated,  with  blue  boxes  denoting  areas  of  glacial  abrasion 
and  green  boxes  denoting  areas  of  glacial  plucking  in  Fig.  5. 

From  the  Larvik,  Norway,  trial,  estimates  of  scattering  strength  from  abraded  surfaces 
range  between  -33  and  -26  dB,  and  scattering  strength  from  quarried  surfaces  range  between 
-30  and  -20  dB.  Scattering  strength  results  are  presented  below,  and  compared  to  model 
curves.  The  measured  scattering  cross  section  from  the  plucked  surfaces  exhibited  variability 
on  the  order  of  5  dB.  The  abraded  part  of  these  outcrops  has  very  low  small-scale  roughness, 
and  it  is  likely  that  the  scattering  cross  section  measured  on  these  sections  can  be  predicted 
by  first-order  scattering  models,  such  as  the  small-slope  approximation  (SSA).  The  plucked 
component  of  these  outcrops  has  very  large  rms  roughness  compared  to  the  wavelength,  and 
it  is  likely  that  first-order  scattering  models  cannot  predict  the  scattering  cross  section. 

From  the  Larvik,  Norway,  trial,  scattering  strength  from  rocks  was  extracted  from  the 
normalized  pressure  squared  by  selecting  a  region  and  averaging  in  cross-range,  and  then 
averaging  over  one  degree  increments.  To  measure  the  scattering  strength  from  a  rock 
surface,  the  mean  slope  was  determined  from  high-resolution  interferometric  bathymetry 
so  that  the  global  grazing  angle  of  the  ideal  mean  seafloor  could  be  mapped  to  the  local 
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Figure  4.  (color  online)  SAS  images  of  the  calibration  rock  outcrop.  Boxes 
denote  areas  where  pixels  were  extracted  to  form  estimates  of  scattering 
strength.  With  respect  to  Ra,  the  rest  of  the  images  (Rb-Rf  in  order)  are 
related  in  azimuth  angle  by  the  following  counterclockwise  rotations:  180°, 
45°,  225°,  270°,  and  390°.  Grayscale  value  denotes  the  decibel  equivalent  of 
a,  the  unaveraged  scattering  cross  section.  Horizontal  axes  represent  ground 
range  from  the  sonar  track  with  positive  values  representing  the  port  side, 
and  negative  values  representing  the  starboard  side.  Vertical  axes  represent 
distance  along  the  sonar  track. 
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Figure  5.  (color  online)  SAS  images  of  roches  moutonees.  Horizontal  axes 
show  gronnd  range  distance  from  the  sonar  track,  with  negative  valnes  cor¬ 
responding  to  port,  and  positive  corresponding  to  starboard.  Vertical  axes 
denote  distance  along  the  sonar  track.  Grayscale  valne  corresponds  to  the 
decibel  eqnivalent  of  d,  the  unaveraged  scattering  cross  section,  and  boxes 
outline  areas  where  pixels  were  used  to  estimate  scattering  strength.  Solid 
boxes  indicate  regions  of  glacial  abrasion  and  dashed  boxes  indicate  regions  of 
glacial  plucking. 
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Figure  6.  (color  online)  Scattering  strength  as  a  function  of  grazing  angle 
from  glacially  abraded  surfaces  a)  R,  the  calibration  rock,  b)  SI,  c)  S2,  d)  S3, 
and  e)  S4.  Error  bars  represent  the  measurement  uncertainty,  and  the  two 
gray  dashed  lines  indicate  the  error  bounds  of  the  small  slope  approximation 
(SSA)  due  to  uncertainty  in  input  parameters. 


grazing  angle  of  the  rock.  Scattering  strength  results  are  presented  in  Fig.  6  for  abraded 
surfaces  and  in  Fig.  7  for  plucked  surfaces.  After  system  calibration,  scattering  strengths 
from  the  low-amplitude  roughness  of  the  abraded  side  of  the  rock  outcrops  was  found  to 
range  between  -33  and  -26  dB  at  20  grazing  angle,  and  agrees  with  predictions  using  the 
SSA  with  input  parameters  measured  during  the  May  2013  experiment.  This  agreement  is 
expected,  since  the  rms  height  of  abraded  surfaces  times  the  acoustic  wavenumber  is  less  than 
unity.  Scattering  strengths  from  the  high-amplitude  roughness  of  the  plucked  side  ranged 
between  -30  and  -20  dB  at  20  grazing  angle.  Predictions  using  small  slope  with  measured 
parameters  from  the  plucked  interface  resulted  in  an  overestimate  of  8  dB  or  more.  This 
overestimate  of  scattering  strength  by  a  first-order  model  is  unexpected,  and  indicates  that 
higher-order  modeling  techniques  may  not  yield  acceptable  predictions  of  scattering  strength, 
due  to  the  fact  that  higher-order  terms  in  the  cross  section  are  typically  positive-definite  [6] . 
The  composite  roughness  model  [7]  is  hypothesized  as  a  plausible  way  forward  to  predict 
the  scattering  cross  section  from  faceted  surfaces  when  using  systems  that  have  resolution 
smaller  than  the  mean  facet  size. 
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Figure  7.  (color  online)  Scattering  strength  as  a  function  of  grazing  angle 
from  glacially  plucked  surfaces  a)  SI,  b)  S2,  c)  S3,  and  d)  S4.  Error  bars  repre¬ 
sent  the  measurement  uncertainty,  and  the  black  line  represents  the  empirical 
Lambertian  model  (LM)  with  the  parameter  p  estimated  by  best-ht  to  the 
data. 


3.3.  Scattering  Statistics  of  rock  outcrops.  The  probability  density  function  (PDF)  of 
the  magnitude  of  the  scattered  held  is  an  important  quantity  for  target  detection  systems. 
This  quantity  is  typically  expressed  in  terms  of  the  probability  of  false  alarm  (PFA),  which 
is  equal  to  1  —  CDF,  where  the  CDF  is  the  cumulative  density  function. 

Analysis  of  the  PFA  from  rock  outcrops  was  performed  during  FY17  and  modeled  using  a 
mixture  model.  The  SAS  images  from  wich  PFA  estimates  were  extracted  are  presented  in 
Fig.  5.  These  images  were  calibrated  during  FY15  and  FY16  of  this  project,  and  scattering 
strength  estimates  were  obtained  and  documented  in  [1]. 

For  this  report,  we  will  restrict  attention  to  the  rock  outcrop  S3  in  Fig.  5.  The  PFA 
is  estimated  by  extracting  pixels  within  the  green  rectangles  in  Fig.  5,  and  binning  them 
according  to  the  large-scale  slope  available  through  SAS  bathymetry.  We  focus  on  angles 
between  20  and  30  degrees  grazing  angle.  The  PFA  from  the  S3  rock  outcrop  is  shown 
in  Fig.  8.  The  data  is  presented  in  the  blue  curve,  and  two  models  are  comapred  to  it,  a 
K-distribution  [8],  and  a  mixture  model  [9].  Although  the  measured  PFA  is  quite  heavy¬ 
tailed,  the  K  distribution  fails  to  capture  the  trends  seen  in  the  data.  The  K-distribution  is 
commonly  used  to  model  patchy  environments  that  behave  like  discrete  scatterers,  and  the 
shape  parameter  is  related  to  the  effective  number  of  scatterers  per  resolution  cell. 

However,  examination  of  S3  in  Fig.  5  shows  that  the  scattered  held  is  caused  by  several 
features  1)  low  amplitude  scattering  from  horizontal  facest,  2)  high  amplitude  from  verti¬ 
cal  facets  that  may  include  multiple  scattering  from  corners,  and  3)  small  dropstones.  If 
we  assume  that  only  the  hrst  two  scattering  sources  dominate  the  scattered  held,  then  a 
two-component  mixture  model  for  the  scattered  held  envelope  distribution  has  a  plausible 
physical  basis.  Measurements  of  the  PFA  of  rock  outcrops  with  low  amplitude  roughness 
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Figure  8.  PEA  of  pixels  bounded  by  green  dashes  of  S3  in  Fig.  5.  The  data 
are  compared  to  both  a  K  distribution  and  a  two-component  mixture  model. 

showed  that  the  low  amplitude  scattering  is  adequately  modeled  using  the  Rayleigh  distribu¬ 
tion.  Since  there  are  a  small  number  of  bright  vertical  facets  in  the  SAS  image,  it  makes  sense 
to  use  a  K-distribution  to  model  the  high-amplitude  component,  because  this  distribution 
results  from  a  small  number  of  discrete  scatterers. 

The  mixture  model  we  use  is  specihed  by 

(1)  p{x\p,  Ao,  a.  A)  =  pPnixlXo)  (1  -  p)pK{x\a,  A) 

where  the  notation  p{y\a,  b,c, . . .)  is  the  pdf  of  random  variable  y  with  model  parameters 
a,b,c, . . . ,  p{x\p,  Ao,  a,  A)  is  the  two-component  mixutre  model  with  four  parameters,  p  is 
the  relative  proportion  of  samples  folloing  Rayleigh  statistics,  Aq  is  the  scale  parameter  of 
the  Rayleigh  distribution,  a  is  the  shape  parameter  of  the  K-distribution,  and  A  is  the 
scale  parameter  of  the  K-distribution.  This  two-component  mixture  distribution  represents 
a  combination  of  background  (Rayleigh)  and  clutter  (K)  distrubitons,  and  was  used  to  model 
midfrequency  active  reverberation  data  in  [9]. 

The  standard  method  to  extimate  mixture  model  parameters  is  called  the  Expectaction 
Maximization  (EM)  algorithm  [9].  The  EM  algorithm  proceeds  in  two  stages.  First  if  we 
have  an  initial  guess  of  the  model  parameters,  {p,  Aq,  d.  A},  we  can  use  those  to  assign  weights 
to  each  data  sample  via 


(2) 

IK  - 

pLo^i  -|-  (1  —  pLi^i 

(3) 

hFi,  =  1  -  hFo, 

(4) 

Lo,i  =  pni^Xp  Ao) 

(5) 

Li,i  =  pK{Xi]a,\) 
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where  hho,*  is  a  weight  that  represents  the  belief  that  sample  i  belongs  to  the  background 
Rayleigh  distribution,  hhi  *  is  a  weight  that  sample  i  belongs  the  clutter  K  distribution,  Lo,i 
is  the  likelihood  function  of  the  background  component  for  sample  i,  Li  i  is  the  likelihood 
function  of  the  clutter  component  for  sample  i,  and  Xi  is  the  ith  data  point.  The  weights 
provide  a  soft  partition  of  the  data  into  each  of  the  mixture  component  distributions.  The 
second  step  is  to  use  these  belief  weights  to  improve  our  estimate  of  the  mixture  parameters 
by  maximizing  an  intermediate  function  called  the  Q  function  for  each  component,  which  is 
the  sum  of  the  product  of  the  log  likelihood  function  for  each  distribution  components  and 
the  weights.  The  Q  function  for  the  background  component  is 

n  n 

(6)  g(p,  Ao)  =  Wo,  log(p)  +  ITo,  log(pn(X„  A)) 

i=l  i=l 

and  the  Q  function  for  the  clutter  component  is 

n 

(7)  g(a.  A)  =  log(pK(X„  a,  A)) 

i=l 

Maximizing  the  Q  function  for  the  background  Rayleigh  component  results  in  a  closed  form 
solution  for  p  and  Aq, 

1  ” 

(8)  p^-TWo,, 

i=l 
1  ” 

(9)  A=— VfRoyX, 

no 

'  2=1 

Maximizing  the  Q  function  for  the  clutter  component  typically  requires  a  numerical  algo¬ 
rithm,  a  few  of  which  are  suggested  in  [9].  In  this  work  for  the  K  clutter  component,  we  use 
the  weighted  method  of  moments  estimator  instead  of  directly  maximizing  the  Q  function  (as 
performed  by  [9]),  and  thus  we  call  the  EM  algorithm  employed  in  this  work  the  EM-MoM. 
This  discrepancy  results  in  different  parameter  estimates  between  the  Maximum  Likelihood 
and  EM-MoM  estimates,  discussed  below. 

These  two  steps  (estimating  the  weights,  and  estimate  the  parameters)  are  repeated  un¬ 
til  some  convergence  criterion  is  met.  The  EM  algorithm  (software  provided  by  Douglas 
Abraham  of  CausaSci  LLC)  was  used  to  estimate  model  parameters  that  are  plotted  in 

Fig.  8. 

The  sequence  of  paremter  estimates  provided  by  EM  can  be  proven  to  monotonically 
increase  the  likelihood  function  [10],  although  such  a  proof  does  not  exist  for  EM-MoM. 
Such  iterative  algorithms  can  often  become  trapped  in  local  maxima,  and  may  not  converge 
to  the  global  maximum  of  the  likelihood  function. 

Apart  from  the  shortcomings  of  the  EM-MoM  algorithm  listed  above,  the  parameter  es¬ 
timates  for  the  proposed  mixture  model  can  be  thought  of  as  random  variables  that  are 
related  to  the  random  variables  that  describe  the  roughness  properties  of  the  rock  outcrops. 
Although  the  roughness  of  the  rock  outcrop  is  deterministic,  it  is  convenient  to  think  of  the 
roughness  properties  as  random.  Since  the  mixture  model  parameters  can  be  linked  via  a 
scattering  model  to  the  roughness  parameters,  it  may  be  possible  to  use  the  scattered  held 
envelope  distribution  as  an  observable  to  remotely  sense  roughness  properties  of  fractured 
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Table  1.  Parameter  estimates  using  expectation  maximization,  and  maxim  a  posteriori 


Parameter 

EM  Estimate 

MAP  Estimate 

P 

0.6526 

0.6066 

-^0 

0.2416 

0.2001 

a 

0.4438 

0.7533 

A 

5.4046 

2.5598 

CBR 

9.9304 

9.6367 

rock  surfaces.  For  use  in  an  inversion  algorithm,  it  is  vital  to  understand  the  joint  distribu¬ 
tion  of  the  roughness  parameters  (which  is  linked  to  the  joint  distribution  of  mixture  model 
parameters)  to  uncover  correlations  and  possible  degeneracies  in  the  data.  Although  EM 
gives  us  parameter  estimates,  it  does  not  tell  us  any  information  about  the  joint  distribution 
of  the  mixture  model  parameters. 

A  Bayesian  framework  can  provide  such  an  estimate  of  the  joint  probability  distribution 
of  mixture  model  parameters.  The  posterior  probability  distribution  (PPD)  of  the  model 
parameters  given  the  data  is  p(m,  d),  where  m  is  a  four  component  vector  of  model  param¬ 
eters,  and  d  is  an  n  component  vector  of  data  points.  The  PPD  is  the  quantity  we  wish  to 
estimate,  and  is  related  to  the  likelihood  function  by 

L(d|m)p(m 
p(m,d)  =  - - - 


where  p(m)  is  the  prior  distribution  of  the  model  parameters  (we  assume  a  bounded  uniform 
distribution),  L(d|m)  is  the  likelihood  function,  and  Z  is  the  evidence,  equal  to  the  prior 
distribution  of  the  data  p(d).  The  evidence  normalizes  the  numerator  in  Eq.  (10)  and  ensures 
that  the  PPD  integrates  to  unity  over  the  parameter  space.  The  best  parameter  estimate  in 
a  Bayesian  framework  is  the  parameter  set  that  maximizes  the  PPD,  and  is  referred  to  as 
the  maximum  a  posteriori  (MAP)  estimate.  If  the  prior  distribution  is  uniform  and  overlaps 
with  the  ML  estimate,  then  the  MAP  and  ML  parameter  estimates  are  the  same.  The  use 
of  bounded  uniform  prior  distributions  for  the  model  parameters  was  motivated  by  both 
our  lack  of  specihc  knowledge  about  the  environment,  and  for  comparing  the  EM  estimates 
with  ML  estimates.  Bounds  of  these  uniform  prior  distributions  have  been  obtained  by 
examination  of  the  data. 

The  PPD  from  the  data  set  in  Fig.  8  was  estimated  using  an  exhaustive  technique,  using 
100  samples  between  the  bounds  of  p  G  {0,0.8},  Aq  G  {0.0120.26},  a  G  {9,1.6},  and 
A  G  {0,6}.  The  exhaustive  search  resulted  in  parameter  estimates  that  can  be  found  in 
Table  1,  and  are  compared  to  the  EM-MoM  estimates.  The  MAP  was  found  by  taking 
the  maximum  over  the  100  sample  4D  parameter  space  as  a  coarse  estimate.  The  coarse 
estimate  was  rehned  using  a  downhill  simplex  algorithm  that  terminated  when  the  log- 
likelihood  function  increased  by  less  than  10“®  from  the  previous  estimate.  The  PPD  in 
this  case  is  a  four-dimensional  distribution  and  is  not  easily  visualized.  Instead,  the  2D 
marginal  distributions  of  the  parameters  are  plotted  in  Fig.  9,  where  each  subfigure  plots 
a  different  combination  of  the  integrated  2D  marginal  distribution.  The  MAP  parameter 
estimates  are  shown  as  green  solid  lines,  and  the  EM  estimates  are  shown  as  dashed  blue 
lines.  Signihcant  correlations  are  seen  between  the  parameters  of  the  mixture  distribution. 
The  one-dimensional  marginal  distributions  are  shown  in  Fig.  10. 
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Figure  9.  Joint  marginal  posterior  probability  distribution  of  the  data  in  S3. 
These  hgures  represent  the  full  PPD  integrated  over  two  parameters,  and  is  a 
convenient  way  to  look  at  the  full  PPD  with  reduced  dimensionality  that  we 
can  visually  interpret. 


The  resulting  model  PFAs  from  these  estimates  are  plotted  graphically  in  Fig.  11.  The 
Rayleigh  component  of  the  mixture  model  is  quite  similar  in  both  the  EM  and  MAP  esti¬ 
mates,  although  the  tails  are  signihcantly  different.  The  different  estimates  in  the  tails  of 
the  distribution  are  conhrmed  by  the  different  estimates  of  the  shape  parameter  for  the  EM- 
MoM  estimate  (about  0.5),  and  the  MAP  estimate  (about  0.75).  In  the  graphs  of  the  PFAs 
in  Fig.  11,  the  mixture  model  tails  seem  to  be  better  estimated  by  EM  than  by  MAP.  The 
“clutter  power”  or  variance  of  the  clutter  component  for  the  K  distribution  is  aA,  and  the 
“background  power”  or  variance  of  the  background  component  is  Aq.  Therefore  the  clutter 
to  background  power  ratio  (CBR)  is  aA/Ao,  and  is  also  included  in  Table.  1.  For  both  EM 
-MoM  and  MAP  estimates,  the  CBR  is  quite  similar,  even  though  the  K  distribution  shape 
parameters  are  different. 

We  hypothesize  that  difference  in  estimates  of  a  and  A  are  due  to  the  method  of  moments 
estimator  used  in  the  maximization  step  of  expectation  maximization.  The  method  of  mo¬ 
ments  is  not  guaranteed  to  result  in  the  same  parameter  estimates  as  maximum  likelihood. 
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Figure  10.  One  dimensional  marginal  distributions  of  the  parameters  of  the 
mixture  model. 


S03 


Figure  11.  PFA  plots  of  the  envelope  data  from  S3  compared  to  mixture 
model  PFAs  with  parameters  estimated  using  the  Expectation  Maximization 
algorithm,  and  MAP  estimate  from  the  Bayes  framework.  Note  that  the 
Rayleigh  components  of  the  both  model  estimates  are  similar,  but  the  tails 
are  signihcantly  different. 


and  the  repeated  use  of  MoM  over  many  iterations  of  the  EM  algorithm  may  cause  the 
EM-MoM  to  diverge  from  the  standard  EM  algorithm. 
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Therefore,  although  the  EM-MoM  algorithm  does  not  necessarily  maximize  the  true  likeli¬ 
hood  function  of  the  mixture  model,  it  provides  a  better  estimate  of  the  tails  and  at  provides 
a  similar  CBR  (at  least  for  this  example).  In  target  detection  systems,  the  distribution  tails 
are  the  most  problematic  and  cause  a  large  number  of  false  alarms.  For  inversion  of  rough¬ 
ness  parameters  using  scattering  statistics,  the  Bayes  framework  may  be  the  most  useful, 
although  it  is  computationally  expensive.  For  modeling  false  alarms  in  active  sonar  systems, 
the  EM-MoM  algorithm  may  be  the  best  method  to  estimate  parameters,  and  it  is  extremely 
fast  compared  to  the  fnll  inversion  for  the  PPD. 

Fntnre  work  in  this  area  will  focus  on  relating  mixtnre  model  parameters  to  roughness 
statistics  of  the  rock  ontcrop,  including  facet  size  distribntion,  and  small-scale  ronghness 
power  spectrum  parameters.  This  physical  model  with  enable  forward  modeling  of  target 
detection  systems,  performance  prediction,  and  potential  for  remote  sensing  of  the  rough¬ 
ness  properties  of  fractnred  rock  ontcrops.  Comparisons  of  the  mixture  model  parameter 
estimates  using  full-resolution  SAS  data  and  reduced  resolution  images  that  simnlate  tar¬ 
get  detection  systems  with  poor  angnlar  resolution  are  key  for  modeling  the  performance  of 
target  detection  systems. 


4.  Impact  and  Applications 

The  primary  work  completed  over  the  course  of  this  project  consisted  of  developing  tech¬ 
niques  for  modeling  scattering  from  rough  rock  outcrop  areas  and  comparing  resnlts  with 
acoustic  data  sets  collected  from  rocky  areas.  The  proposed  project  was  designed  to  increase 
our  understanding  of  and  simulation  capability  for  acoustic  scattering  from  rock  outcrops. 
This  study  yielded  useful  knowledge  of  rock  ontcrops  as  a  mechanism  responsible  for  shallow 
water  false  alarms,  and  how  levels  of  false  alarm  relate  to  physical  properties  and  features 
of  rock  outcrops.  Guidance  relevant  to  this  false  alarm  mechanism  is  being  provided  to 
researchers  at  the  Naval  Research  Laboratory  and  will  also  be  provided  to  those  developing 
digital  simulation  content  as  part  of  other  programs.  Other  deliverables  are  journal  articles 
that  have  been  published,  and  are  in  preparation. 

We  collaborated  with  researchers  at  NRL  on  their  related  project,  “Modeling  of  High- 
Freqnency  Broadband  False  Target  Phenomen”  (project  Pis  were  Roger  Gauss,  Dave  Galvo, 
and  Joe  Fialkowski)  as  well  as  continning  the  joint  research  project  “Gharacterization  and 
Modeling  of  Synthetic  Aperture  Sonar  Data”  with  FFI  (POG  -  Roy  Hansen).  The  collabo¬ 
ration  with  NRL-DG  resulted  in  two  conference  publications  that  will  be  presented  at  the 
OGEAN15  conference  in  October  2015  [11,  12]. 
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